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Abstract 

The dependence between neutron-to-proton yield ratio (R np ) and neutron skin thickness (6 np ) in neutron- 
rich projectile induced reactions is investigated within the framework of the Isospin-Dependent Quantum 
Molecular Dynamics (IQMD) model. The density distribution of the Droplet model is embedded in the 
initialization of the neutron and proton densities in the present IQMD model. By adjusting the diffuseness 
parameter of neutron density in the Droplet model for the projectile, the relationship between the neutron 
skin thickness and the corresponding R np in the collisions is obtained. The results show strong linear 
correlation between R np and S np for neutron-rich Ca and Ni isotopes. It is suggested that R np may be used 
as an experimental observable to extract 6 np for neutron-rich nuclei, which is very significant to the study 
of the nuclear structure of exotic nuclei and the equation of state (EOS) of asymmetric nuclear matter. 
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Nuclear radius is one of the basic quantities of a 
nucleus. The proton root-mean-square (RMS) ra- 
dius can be determined to very high accuracy via 
the charge radius measured by electromagnetic inter- 
actions, typically with an error of 0.02 fm or better 
for many nuclei In contrast, it is much more 
difficult to accurately determine the neutron density 
distribution of a nucleus experimentally Q. Thus, 
the accuracy of experimental neutron radius is much 
lower than that of the proton radius. However, the 
information of neutron density is very important to 
the study of nuclear structure for neutron-rich nu- 
clei, atomic parity non-conservation, iso- vector inter- 
actions, and neutron-rich matter in astrophysics etc. 
It is remarkable that a single measurement has so 
many applications in the research fields of atomic, 
nuclear and astrophysics d, 0] . 
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A nucleus with neutron number (N) being larger 
than proton number (Z) are expected to have a neu- 
tron skin (defined as the difference between the neu- 
tron and proton RMS radii: 8 np = (r 2 ) 1 / 2 - (r 2 ) 1/2 ). 
The neutron skin thickness 5 np depends on the bal- 
ance between various aspects of the nuclear force. 
The actual proton and neutron density distributions 
are determined by the balance between the isospin 
asymmetry and Coulomb force. 6 np is found to be 
related with a constraint on the equation of state 
(EOS) of asymmetric nuclear matter. Strong linear 
correlation between 5 np and L (the slope of symme- 
try energy coefficient C sym ), the ratio L/J (J is the 
symmetry energy coefficient at the saturation density 
p ), J - a sym (a sym is the symmetry energy coeffi- 
cient of finite nuclei) are demonstrated [B[ . This con- 
straint is important for extrapolation of the EOS to 
high density and hence useful for studying properties 
of neutron star @, 0, 0, S, H, E3, EH, E3 ■ Neutron skin 
thickness can yield a lot of information on the deriva- 
tion of volume and surface symmetry energy, as well 
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as nuclear incompressibility with respective to den- 
sity. S np is significant to the study of the EOS in 
different theoretical models such as Skyrme Hartree- 
Fock (SHF) [| SJ3J1 & EH, El, relativistic mean- 
field (RMF) EJmlll, BUU model 0, EE El and 
Droplet model |l7| . Furthermore, neutron skin thick- 
ness helps to identify a nucleus with exotic structure. 
Thus the precise determination of S np for a nucleus 
becomes an important research subject in nuclear 
physics. 

Several attempts have been made or suggested to 
determine the neutron density distribution such as 
using proton scattering, interaction cross sections in 
heavy ion collisions at relativistic energies [l8| , parity 
violating measurements 0, 0, 0, S, Q , neutron abra- 
sion cross sections in heavy ion collisions On 
the other hand, the neutron and proton transverse 
emission and double neutron-proton ratios have been 
studied as a sensitive observable of the asymmetry 
term of the nuclear EOS in the experiment [2011 and 
different kinds of simulations 2l|, |22j, |23j, ]2M, \2M ■ ^ n 
this Letter, the relationship between 5 np and the ra- 
tio of the emitted neutron and proton yields (R np = 
Y n /Yp) in neutron-rich projectile induced reactions 
is, for the first time, presented within the framework 
of Isospin-Dependent Quantum Molecular Dynamics 
(IQMD) model. The possibility of extracting S np 
from R np is investigated. 

The QMD approach is a many-body theory that 
describes heavy-ion reactions from intermediate to 
relativistic energies [26j. It includes several impor- 
tant factors: initialization of the projectile and tar- 
get, nucleon propagation in the effective potential, 
nucleon-nucleon (NN) collisions in a nuclear medium 
and the Pauli blocking effect. A g eneral review of the 
QMD model can be found in [27] . The IQMD model 
is based on the QMD model with consideration of the 
isospin effect. 

The dynamics of heavy ion collision at intermediate 
energies is governed mainly by three components: the 
mean field, two-body collisions, and Pauli blocking. 
Therefore, for an isospin-dependent reaction dynam- 
ics model, it is important to include isospin degrees 
of freedom with the above three components. In ad- 
dition, the sampling of phase space of neutrons and 
protons in the initialization should be treated sepa- 



rately because of the large difference between neu- 
tron and proton density distributions for nuclei that 
far from the /3-stability line. 

In the IQMD model, each nucleon i is represented 
by a Gaussian wave packet with definite width (L — 
2.16 fm 2 ) centered around the mean position and the 
mean momentum: 
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The nuclear mean field in the IQMD model is deter- 
mined as following 
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with the normal nuclear matter density po = 
0.16 fm~ 3 . p, p n and p p are the total, neutron and 
proton densities, respectively. t z is z-th component 
of the isospin degree of freedom, which equals 1 or 
— 1 for neutron or proton, respectively. The coef- 
ficients a, [3 and 7 are parameters for the nuclear 
EOS. C sym is the symmetry energy strength due to 
the difference between neutron and proton. In the 
present work, we take a = —356 MeV, f3 = 303 MeV 
and 7 = 1.17 which correspond to the so-called soft 
EOS with C sym = 32 MeV and incompressibility of 
K = 200 MeV. V c is the Coulomb potential and U Ynk 
is Yukawa (surface) potential. In the phase space ini- 
tialization of the projectile and target, the density 
distributions of proton and neutron are distinguished 
from each other. The neutron and proton density 
distributions for the initial projectile and target nu- 
clei in the present IQMD model are taken from the 
Droplet Model 

In the Droplet model, we can change the diffuseness 
parameter to get different density, 
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where p® is the normalization constant which ensures 
that the integration of the density distribution equals 
to the number of neutrons (i=n) or protons (i=p); U 
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Figure 1: (Color online.) The evolution time dependence of 
Rnv for 50 Ca+ 12 C at 50A MeV. 



is the diffuseness parameter; Ci is the half density ra- 
dius of neutron or proton determined by the Droplet 
model 12811 



C i = R i [l-{b i /R i )%i = n,p 



(4) 



here hi — 0.413/^, Ri is the equivalent sharp surface 
radius of neutron and proton. Ri and ti are given by 
the Droplet model. A factor /j is introduced by us to 
adjust the diffuseness parameter. In [29], Trzci'nska 
et al. found that the half density radii for neutrons 
and protons in heavy nuclei are almost the same, but 
the diffuseness parameter for neutron is larger than 
that for the proton. So we introduce the factor /j to 
adjust the diffuseness parameter for neutron. In the 
calculation for neutron-rich nucleus, f p — 1.0 is used 
in Eq.(3) for the proton density distribution as in the 
Droplet model, while f n in Eq.(3) is changed from 1.0 
to 1.6. Different values of S np will be deduced from 
Eq.(3). Using the density distributions of the Droplet 
model, we can get the initial coordinate of nucleons in 
nuclei in terms of the Monte Carlo sampling method. 
The momentum distribution of nucleons is generated 
by means of the local Fermi gas approximation: 
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In the IQMD model, the nucleon's radial density can 
be written as 
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Figure 2: (Color online.) The dependence of R np on neu- 
tron skin thickness of the projectile for 50 Ca+ 12 C (a) and 
6s Ni+ 12 C (b) at 50A MeV under the condition of Y > 0. 
Different symbols are used for different range of the reduced 
impact parameter as shown in the legend. The dotted lines 
just guide the eye. 



To avoid taking an unstable initialization of pro- 
jectile and target in the IQMD calculation, we only 
select the initialization samples of those nuclei that 
meet the required stability conditions. The bind- 
ing energies and root mean square (rms) radii of 
the initial configuration for the colliding nuclei are 
required to be stable. Using the selected initializa- 
tion phase space of nuclei in IQMD to simulate the 
collisions, the nuclear fragments are constructed by 
a modified isospin-dependent coalescence model, in 
which nucleons with relative momentum smaller than 
Pq = 300 MeV/c and relative distance smaller than 
Rq = 3.5 fm will be combined into a cluster. 

The collision processes of some Ca and Ni isotopes 
with 12 C target at 50 A MeV are simulated using the 
IQMD model. The fragments including neutrons and 
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Figure 3: (Color online.) The dependence of R np on neutron 
skin thickness under the condition of reduced impact param- 
eter from 0.8 to 1.0 for 50 Ca+ 12 C (a) and 68 Ni+ 12 C (b) at 
50A MeV. The results without any gate on the rapidity (Y) 
are shown as solid squares. The results with Y > are plotted 
as solid circles. The dotted lines just guide the eye. 



protons that formed during the evolution of the colli- 
sion are constructed by the coalescence method. The 
yield ratio R np of the emitted neutrons and protons 
can be calculated from the yields of the produced 
neutrons and protons. By changing the factor /„ in 
the neutron density distribution of the Droplet model 
for the projectile, different values of 5 np and the cor- 
responding R np are obtained. Thus we can obtain 
the correlation between R np and S np . In the calcula- 
tion, the time evolution of the dynamical process was 
simulated until t = 200 fm/c. The calculated R np is 
stable after 150 fm/c as shown in Fig. [TJ so we ac- 
cumulate the emitted neutrons and protons between 
150 fm/c and 200 fm/c in order to improve statis- 
tics. The R np from different range of the reduced 
impact parameter for 50 Ca + 12 C and 68 Ni + 12 C are 



plotted in Fig. [3] The reduced impact parameter is 
defined as b/b max with 6 max being the maximum im- 
pact parameter. From Fig. [21 we can see that R np 
rises as S np increases. With S np being fixed, R np also 
becomes larger with the increasing of the reduced im- 
pact parameter. It means that R np from peripheral 
collisions is usually larger than that from central col- 
lisions. The main purpose of the present study is 
to investigate the relationship between R np and the 
neutron skin thickness of the projectile. In order to 
minimize the target effect on R np , we use rapidity (Y) 
cut to select neutrons and protons from the projec- 
tile. The rapidity of the fragment which is normalized 
to the incident projectile rapidity is defined as: 

where E is the energy of the fragment, p z is the mo- 
mentum of z direction, Y pm j is the rapidity of the 
projectile. We choose Y > to strip away fragments 
coming from the target. From the results shown in 
Fig.[3J we can see that R np with Y > are larger than 
R np without Y cut. Since the projectile is neutron- 
rich nucleus, the correlation between R np and 6 np is 
stronger for nucleons coming from the projectile than 
those from both the projectile and target. 

Strong linear correlation between R np and 5 np of 
the projectile, especially for 0.8 < 6/6 max < 1.0, is 
exhibited in Fig. [2] It indicates that R np is very 
sensitive to 5 np of the projectile, especially in pe- 
ripheral collisions. For systematic study, reactions 
of other Ca and Ni isotopes such as 52 > 54 > 56 Ca and 
60,62,64,66,70 Ni are als0 calculated. The results with 
the reduced impact parameter from 0.8 to 1.0 and the 
rapidity being positive are shown in Fig. [4] Strong 
linear relationship between R np and S np are also ob- 
served for all projectiles. From Fig. [4j a linear func- 
tion (R np — a + b ■ S np ) can describe the correlation 
between R np and S np well for both Ca and Ni iso- 
topes. The dependence between R np and 5 np is fit- 
ted using this linear function. The mean slopes for 
Ca and Ni isotopes are 1.06 and 0.83, respectively. 

From the above discussions, R np could be viewed 
as a sensitive observable of S np for the projectile. If 
the produced neutrons and protons are measured by 
experiment, it is possible to extract S np from the 
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Figure 4: (Color online.) The dependence of R np on neutron 
skin thickness under the condition of Y >0 and 0.8 < b/b ma , x < 
1.0 for 50,52,54,56 Ca+ 12 c (a) and 60,62,64,66,70 Ni+ 12 c (b) at 

50A MeV. Different symbols are used for different projectile as 
shown in the legend. The dotted lines just guide the eye, for 
detail see text. 
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Figure 5: (Color online.) The dependence of R np on the 
mass number A under the condition of Y > and 0.8 < 
&/6max < 1.0 for 50,52,54,56 Ca+ 12 G and 60,62,64,66,70 Ni+ 12 c 

at 50A MeV. Different symbols are used for different value of 
the factor /„ as shown in the legend. The dotted lines just 
guide the eye, for detail see text. 



ncutron-to-proton yield ratio R np . From the fitted 
slope values, the estimated error of S np may be less 
than 0.08 fm if the uncertainty of R np is less than 5%. 
It should be pointed out that this uncertainty does 
not take into account the error arising from the de- 
termination of impact parameter, which is always an 
important issue in studies of nuclear reactions in the 
energy range of our calculation. Wc estimate the un- 
certainty of S np for 0.6 < 6/6 max < 1.0 to be around 
0.1 fm when the uncertainty of R np is 5%. To achieve 
high resolution, well defined experimental determina- 
tion of the centrality of the reaction is required for 
this method. 

In our calculation, different S np is obtained by 
changing the factor /„ in Eq.(3). We can also study 
the mass dependence of R np with the same f n as 
shown in Fig. For /„ = 1.0, it refers to the neu- 



tron skin thickness of the Droplet model's prediction. 
In this figure, the linear dependence of R np with the 
mass number A is also seen at different value of /„. 
Due to the acceptance and efficiency in measurement 
of the emitted neutron and proton, the uncertainty 
of R np may be not easy to be estimated which makes 
it difficult to obtain the error of 5 np . There is one 
method to minimize these effects: we can measure 
R np of one isotope chain with different mass number. 
In this case, the acceptance and efficiency of neutron 
and proton will not change the mass dependence of 
R np . If 5 np of some nuclei is larger than the normal 
value of the Droplet model, it will deviate from the 
line of /„ = 1.0. By studying the mass dependence of 
R np , it is possible to extract S np and identify nucleus 
with abnormal neutron skin thickness. 

The neutron excess (/ = (N — Z) /A) dependence of 
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neutron skin thickness for four parameter sets (NL3, 
NL-SH, SkM*, SIII) in the Droplet model and self- 
consistent extended Thomas-Fermi (ETF) calcula- 
tions have been investigated by Warda et al. (30j . 
Differences of the neutron skin thickness among dif- 
ferent predictions increase with the increase of I. To 
distinguish which potential parameter set or method 
is more close to the experimental data, higher reso- 
lution measurement on neutron skin thickness is re- 
quired for stable nuclei having small / in comparison 
with nuclei far from the stability line having large I 
value. From their results, the neutron skin thickness 
by different calculation varies from 0.01 fm to 0.11 
fm at I = 0.11. In this sense, it will be capable of 
distinguishing different theoretical predictions when 
/ > 0.11 by extracting S np from R np measurement 
which the present paper proposes. In other words, 
the present method will be more feasible for heavy 
nuclei or intermediate mass nuclei far from the j3- 
stability line. 

In summary, we have made calculation on the rela- 
tionship between R np and the neutron skin thickness 
for the first time. The simulated data for Ca and Ni 
isotopes induced reactions show a good linear corre- 
lation between R np and S np for neutron-rich projec- 
tile. It is suggested that R np could be used as an 
experimental observable to extract the neutron skin 
thickness for neutron-rich nucleus. When the neutron 
skin thickness of one isotope is obtained, we can get 
some information on the EOS in different theoreti- 
cal models such as Skyrme Hartree-Fock, relativistic 
mean-field and BUU model. From the isotope depen- 
dence of neutron skin thickness, a lot of information 
on the derivation of volume and surface symmetry 
energy, as well as nuclear incompressibility with re- 
spective to density can be extracted. This is very 
important to the study of the EOS of asymmetric 
nuclear matter. 
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